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ABSTRACT 


Conventional approaches for a distributed low probabil- 
ity of detection communications system with a large number 


of unique transmitters, and a single or a few receivers, re- 








quire receiver complexity proportional to the number of 


transmitters. To improve efficiency in terms of receiver 





complexity, two alternative designs are analyzed and com- 








pared to a reference receiver whose complexity grows line- 


arly as the number of transmitters increases. The first al- 





ternative system groups the transmitters into clusters whose 
pseudorandom noise codes have some chips in common. The re- 
sulting receiver would then perform two stages of processes: 


identification of the transmitting cluster and received bit 








detection. The total number of processes required for any 


given transmitter would be substantially less than the tra- 





ditional receiver. The second alternative design would 


utilize a common long spreading code and a shorter cycli- 





cally shifted spreading code in each transmitter. The re- 





ceiver utilizes the cyclic shift property of the fast Fou- 
rier transform to recover efficiently both the identity of 


active receivers and the data sent using a single branch. 





The complexity of the two proposed systems is compared to 


that of the reference system. 
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EXECUTIVE SUMMARY 


Large distributed networks often have more transmitters 





than receivers, like the uplink in a single cellular network 








cell. Such a system can also use direct sequence spread 

















spectrum (DSSS) in order to achieve a low probability of de- 











tection (LPD) and multiple access. DSSS systems achieve LPD 
performance using pseudo-random (PN) codes. Such coding al- 


so improves anti jamming (AJ) performance and mitigates in- 





terference between transmitters. 


All designs considered in this thesis have transmitters 








that use a unique PN code for LPD, AJ, and multiple access. 








In order for a receiver to process all signals, it must have 











memory to stor th referenc PN code(s), a mixer to de- 





spread the signal, a matched filter to extract decision da- 


ta, a sampler to select discrete data for analysis, and a 





decision algorithm in hardware or software to make bit deci- 











sions. In a traditional receiver, this requirement results 
in significant growth in complexity as the number of trans- 
mitters increases, making large systems unwieldy to field 


and difficult to scale for increased traffic loads. 





One potential solution involves sub-dividing the re- 





ceived signal set into smaller groups using a mask of common 





PN code points. For example, a system with 32 transmitters 





requires 32 copies of the receiver processing hardware, or a 


software radio eguivalent, including a demodulation algo- 





rithm run 32 times. Dividing the set into four subsets of 





eight related signals improves the efficiency of processing, 





as a single signal requires only a portion of the receivers 








for demodulation. The receiver would first determine which 


xX1il 


subset contains the received signal, then attempt to demodu- 





late all signals in that subset. The total number of proc- 








esses would be 12 vice 32-four initial detection decisions 








to determine the subset and eight demodulations. This is 





potentially efficient if the number of active transmitters 





is a small fraction of the total transmitters at any given 





moment. However, in the case of all transmitters sending 
nearly simultaneously, the coding mask receiver is less ef- 


ficient as 32 demodulations and four detections require 36 





processes, while the linear-growth receiver still only re- 





quires 32. Additionally, the probability of error in the 





coding mask system will be worse than for the reference r 
ceiver as the mask detection stage introduces additional er- 


ror. 


A second proposed solution involves using short PN 





codes related to a single code by cyclic shift for multiple 


access. To provide for code division multiple access 





(CDMA), each transmitter is preprogrammed with a predeter- 
mined shifted form of a base PN code with duration of a sin- 


gle bit. Every transmitter then uses a common longer PN 





code to provide transmission security by mixing the two 
codes together prior to mixing with data. The longer PN 


code has duration greater than a single bit and potentially 








never repeats within the duration of th ntire message th 





reby ensuring a unique PN code for each transmitter and LPD. 


This proposed receiver also removes the common PN code 





through mixing before filtering and sampling of the combined 





received signal. It compares the discrete Fourier transform 








(DFT) processed signal received against a reference DFT of 








the short unshifted PN code, extracts the transmitter iden- 


Xi1v 











tity through an inverse DFT (IDFT) and completes bit deci- 





sions. This process achieves the same functionality as the 


matched filter in the reference receiver, but with less com- 





ponents required. 








The IDFT output identifies which transmitters sent 





data, as well what data was sent. Increasing the number of 





transmitters requires additional processing of the IDFT out- 














put, but this growth is smaller as the system needs fewer 
copies of shift size and bit detection components. Finally, 


the cyclic PN code system’s improved scalability is not at 








the expense of error performance. It has equivalent error 


performance to the reference receiver for binary phase shift 








keying (BPSK) modulation within an additive white Gaussian 





noise (AWGN) channel. 
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I. INTRODUCTION 


A. BACKGROUND 


This thesis examines potential solutions to the problem 





of scalability for a widely distributed direct sequence 








spread spectrum (DSSS) communications system. This system 





must have low probability of detection (LPD) characteris-— 


tics, which prevent the exploitation of its signals by an 





interceptor, as well as prevent interference between trans- 








mitters. Each transmitter, therefore, must have a unique 





pseudo-random noise (PN) code, and that in turn requires th 
receiver to process all possible system PN codes simultane- 


ously. 


Adding extra copies of the processing hardware for each 


additional transmitter quickly leads to an unwieldy system. 








There will be growing hardware costs and physical space re- 
quirements, as well as potential increased delays in the 
signal processing. Throughout this thesis, the system that 
requires an additional set of hardware for each additional 


transmitter will be referred to as the “reference receiver” 























and its detailed analysis is contained in Chapter : The 





reference receiver will provide the error performance bench- 








mark for comparison of scalability solutions. The referenc 





receiver and all solutions to the scalability problem will 








use DSSS with unique PN codes in each transmitter to main- 





tain LPD and multiple access. Additionally, each analysis 


will use an additive white Gaussian noise (AWGN) channel 





with no forward error correction coding. 





Low probability of interception (LPI) and low probabil- 








ity of exploitation (LPE) are sometimes used interchangeably 





with LPD. LPD in the context of this thesis refers to 











structuring the signal such that it is relatively immune 
from being detected [1], and this term will be used from now 


on in discussing preventing the exploitation of the signal. 


B. OBJECTIVE 





This thesis will examine and validate the reference r 


ceiver, as well as analyze two potential solutions to the 





scalability problem by comparing their error performance and 














scalability. One solution is the “cyclic PN code receiver.” 








It uses two PN codes to provide multiple access and identi- 





fication of transmitters, as well as LPD. A common receiver 


with a single branch of hardware processes all signals si- 








multaneously using discrete Fourier transforms (DFTs), which 





act as a matched filter with a single bit decision block. 


An alternative solution is the “coding mask receiver.” 








It uses the same hardware as the reference receiver for bit 





decisions but sub-divides them into groups with similar PN 
codes. Each subset then has a detection branch to reduce 


the total number of processes required for a single received 





signal. This solution speeds up processing of communica- 


tions by requiring fewer processes for a single signal re- 





ceived. However, the same physical hardware is required to 





support the full system, so scalability is still potentially 


an issue. 


Cc. RELATED WORK 


PN codes with different lengths in conjunction with the 





discrete Fourier transform provide a method to efficiently 
emulate a set of matched filters and maintain bit error rate 


performance. The cyclic PN code receiver solution is simi- 





lar to cyclic code shift keying (CCSK), a technique used in 














the Joint Tactical Information Distribution System (JTIDS) 





[2], except that CCSK uses the cyclic shift to indicate a 


data value while in this thesis, the cyclic shift provides 








multiple access and identifies the source transmitter at the 


receiver [3]. The coding mask receiver solution uses unique 





PN codes with points of similarity to perform the initial 





signal detection and then separate subgroups of receivers 


for demodulation. In terms of the signals, this is the same 








as a typical CDMA [4], but with prescribed correlation be- 





tween portions of the spreading PN codes. 





Additionally, the cyclic PN code receiver’s shorter PN 


code repeats at the bit rate and provides code division mul- 





tiple access (CDMA) similar to the methods described by Lee 





and Miller and separately by Ha as used in the IS-95B cellu- 








lar standard. There, the long PN code scrambles the data 


while two short codes, with different shifts, maintain mul- 








tiple-access within and between cells while minimizing in- 








terference [4], [5]. 


The next chapter provides an overview of the design 


points of the various receivers investigated in this thesis. 




















Chapter then provides a detailed analysis and perform—- 








ance simulation of the reference receiver. Chapter IV fol- 








lows with a detailed description, theoretical calculations, 


and simulation of the cyclic-PN code receiver, followed in 
3 


Chapter V by a Similar analysis of the coding mask receiver. 








Chapter VI provides a summary of all data and compares error 


rate and scalability performance. Finally, conclusions are 





provided in Chapter V 











II. SYSTEM DESIGN 











In order to operate as a distributed communications 
system, a proposed design must meet thr requirements. TSG 
must maintain low probability of detection (LPD). This 








means preventing the detection of the signal by an intercep- 








tor through the design of the signal’s waveform and does not 
refer to protecting the data through encryption. Next, a 
proposed design must have error performance similar to un- 
coded BPSK in an AWGN channel. Superior scalability 
achieved through an increase in the bit error rate is not 


acceptable. Finally, each valid solution must have superior 





scalability to the reference receiver, which is expected to 


have linear growth as the number of transmitters increases. 





A. LOW PROBABILITY OF EXPLOITATION DESIGN 











Direct sequence spread spectrum (DSSS) addresses issues 





of LPD by spreading a signal’s spectrum to decreas th 





probability of detection [1]. Additionally, using multiple 
pseudo-random noise (PN) codes allows for increased protec-— 
tion as well as improvements in demodulation as each re- 


ceiver can use unique inner codes. 


1. Non-recursive Coding 


Effective PN codes must have a random appearance. Each 


transmitter must therefore us a non-recursiv code, 1.e., 








one that does not repeat for the duration of each transmis-— 





sion and ideally never use the same code in subsequent 
transmissions. Mixing two PN codes achieves this as long as 


one of the codes is non-recursive [8]. 


5 


A random number generator, as simulated in MATLAB, can 





generate an infinite length code for this purpose. Ina si- 


mulation, a common code generator can provide codes for the 





transmitter and receiver. In a real system, synchronizing 





the generators with a common seed achieves the same [9], 


[10]. 


2. Low Power 


An LPE communications system typically uses less power 


in its transmitters in order to reduce the chance of detec- 














tion. However, without forward error correction coding, 
this will also reduce error performance and range. For the 
models in this thesis, the energy per bit remains constant 


while the noise power varies, allowing for the testing of a 


range of signal to noise ratios. 


3. High Bandwidth 


Spread spectrum signals require more bandwidth than 


non-spread signals to transmit at the same data rate [5]. 





The proposed solutions control mutual interference through 
PN codes. Future design work should explore forward error 
control coding to support binary phase shift (BPSK) or quad- 


rature phase shift keying (QPSK). This thesis’ simulations 








used BPSK modulation for simplicity. However, DSSS designs 


using QPSK have superior LPE performance, since the signal 





squared does not reveal the carrier frequency [1]. 


B. SCALABILITY 


Examining a system with N transmitters, for a linear- 


growth receiver to processes the received signals, it must 





repeat the filtering and bit decisions for all signals WN 


times, resulting in an approximate total number of multipli- 








i 2 : 
cation processes of N°’. However, the fast Fourier trans-— 





form (FFT) algorithm, a means to implement the discrete Fou- 





rier Transform (DFT), requires 0.5Nlog,N complex multiplica- 


tion operations [6]. The proposed receiver will have one 





FFT and one inverse FFT (IFFT) each, requiring a potential 





total of Nlog,N complex multiplications, which increases at 


a rate less than the reference receiver as the system grows 





in size. A detailed scalability discussion is included 








Chapter V, section B. The next chapter describes the refer 


nce receiver in detail. 
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III. REFERENCE RECEIVER 





The reference receiver uses a unique PN code for each 





transmitter. Its receiver has some common components that 





process all received signals together, but it relies on sep- 


arate de-spreading, filtering, and bit decision blocks for 





each transmitter. The reference receiver is representative 


of a straightforward design approach, and therefore provides 








a good benchmark for comparisons. Determining the referenc 











receiver’s BPSK error performance and validating theoretical 
calculations with simulation are the basis for comparison 


and analysis of the two alternative designs. The block dia- 





gram for this reference design was based on Nicholson’s hy- 








pothetical DSSS transmitter with a single PN code and car- 


rier frequency oscillator [1]. 


A. REFERENCE RECEIVER DESIGN 


Th referenc transmitter and receiver use a single 








spreading code and upconversion to create the spread spec-— 


trum waveform. The receiver complexity grows linearly as 





the number of transmitters increases, but it provides the 





benchmark for error performance. 


2E, 








Its probability of bit error is of } where E, is 


0 





energy per bit and Noy is the two-sided noise power spec 





tral density, which is identical to a binary phase shift 


keyed (BPSK) system [7] and derived in equations (3.16) 


and (3.17). This uncoded BPSK error performance iS aS ex-— 
pected for a DSSS system because the signaling is antipodal 


[8]. 





1. The Reference Transmitter 


The reference transmitter initially transforms the po- 








lar binary data stream {d,\ into am antipodal non-—return-to 


zero (NRZ) binary signal m, (t) at the bit rate R= Yh ; 
b 


«th 


where Pp, (t-iT, ) is a pulse function for the © bit time with 


curation i, and: subscript: yp an m, (t) is the transmitter 


number as shown in equation (3.1). 


m, (t)=>\d,pp,(t-iT,) iG.) 
The data mixes with the transmitter’s unique pseudoran-— 


dom noise (PN) code, CAL), and carrier wave, cos(27f.t) . 








Each transmitter uses the same carrier frequency (f,), which 


is much higher than the bit rate. The amplifier applies a 





preset gain of A before the signal enters the channel as 


s(t)=Am, (t)c, (t)cos(27f.t) , as shown in Figure 1. 


cos (27 f.t) 
m, (t)€ {1,-l} ——>(X)—>(X) > s(t) 


Figure 1. p” Reference Transmitter 
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2. The Reference Receiver 


The channel for this model is assumed to be additive 


white Gaussian noise (AWGN), adding a white Gaussian noise 


factor n(t) to the combined signal, resulting in a received 


signal of s(t)+n(t) as shown in equation (3.2). At the out- 


put of the LNA, the signal received still matches equation 


(3.2) for theoretical calculations. 





r(t) = Am, (t)c, (t)cos(27f,t)+n(t) (3.2) 





One complete receiver branch is required for each of N 


transmitters, resulting in a system as shown in Figure 2. 


cos (27 ft) c[n] 
















Decision > 1 bit 



































LA. Decision -> 1 bit 























Decision > 1 bit 





Figure 2. Reference Receiver System 





Examining just the signal received from the p" trans-— 
mitter, the received signal passes through a low-noise am- 
plifier (LNA) before mixing with the carrier and unique PN 


code matching its transmitter. The signal now has two com-— 





ponents, the squared cosine function, and noise. 
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r(t)=Am, (t)c, (t)c, (t)cos(27f,t)cos (27 f.t)+n(t)c, (t)cos(27f,t) foe 


=Am,(t)c,’(t)cos* (27 f.t)+n(t)c, (t)cos(27 ft) 
Using the trigonometric identity cos" (x) = (I-cos(2x)) 


[11], and the squared PN code equating to one, equation 





(3.3) simplifies, giving a signal with three components: a 


signal term, a double frequency term, and a noise term. 


a(t)=—m, (t)(1-cos(4zf,t))+n(t)c, (t)cos(2z ft) 
(3.4) 


=<m,(t)-Sm, (t)eos(4z ft) +n()¢, (1)e0s(2zf,) 


Looking at only the signal term entering the matched 


filter, its output is 


y(t)=n(r)+4m,(1)=4 | pp (a)m,(t-a)da. (3.5) 


With substitution, the time interval of the convolution 








integral changes from -o<a<o to O<a<iT,. Using the new 





interval, the pulse function term, Pr(a), in the convolution 








integral equates to one, giving a new form of the integral. 


ice} T, 
= 23.4, tp, (t-a-iT, a (3.6) 
oO 0 





Evaluating this integral at the sample instant, t=AT,, 
where keZ, yields the signal portion of the decision sta- 
tistic. 


ioe) 


T, 
kt, )=5 4, pr ((k-i)T, -a)da (3.7) 
—o 0 








In evaluating the convolution integral with k,ieZ, 


there are three possible relations: k-i22, k-i<0, and 








k-i=1. In. the first. imstance, the integral is. zero in 


0<a<iT,. The sécond condition likewise solves to .zéero- for 
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the same reason. In the third instance, evaluating the in- 





tegral at k-i=1 results in the final form of the signal of 





interest at the output of the matched filter. The final 
value shown in equation (3.8) is the mean for the decision 


statistic normal random variable after sampling and prior to 





the bit decision. 
A foo} T, 
y(KT,) =F de, | pr (T,-@)da 
ara - 


=<d, ,| p;,(T,-a@)da (3.8) 
0 


3. Double Frequency Term 


A 
The double frequency term, Me (#)cos(4zf.t) convolves 





with the matched filter impulse response, h(t)=p,(t), to 


yield 


(39) 
=— | p, (@)n(t-a)cos(4z f, (t-a)) da. 


Using the same substitution for the data signal as in 
the signal of interest above, the interval changes from 


—0<a<0o to O<a<T,, giving a new convolution integral of 


the form: 


o TF, 
y()= 4 p, (t-a-iT)cos(4zf,(t-a))da . (3210) 
- 9 


Evaluated at t=kT,, where keZ, the integral becomes: 


ice) 


y(IT,)= 2904, pp (kA) 1, coos arf, (kt, -a))da. (3.11) 


0 
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Similar to the signal of interest, the double frequency 


convolution integral only produces non-zero output when 





k-i=1, resulting in the final form shown in equation 


(3.12). 


Tt, 
y(KT,) =F de] Pr (T,-—a)cos(4zf, (kT, -a))da 
=F 4.1] 00s (4zf, (kT, -a))de (3.12) 


0 


= See sin 4&7) sin ef ((k -1)T,))| 


a c. 
The term in square brackets in the solution has a max 


value of +2 and minimum of -2. Therefore, the absolute 





value of the solution must be: 


 2Adns 


(3:2:13:) 
82 f. 


aan) (4-1) 





Accounting for f,>>l and its location in the denomina- 





tor, all possible results of this term are small compared to 





the possible results of the data signal in equation (3.8), 





so the double frequency term is ignored as insignificant in 


further analysis of the bit decision statistic. 
4. Noise Calculations 


The noise term enters the match filter as 


n(t)c,(t)cos(27f.t) . However, analyzing the resulting power of 





the noise is more useful. The power spectral density (PSD) 
No 

a 
Mixing with the PN code has no effect on the resulting PSD. 





of noise from the AWGN channel at the antenna is Su (t)= 





. . . . 2 
The instantaneous power of a cosine function is cos’ (x), 
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1 
which has an averaged of oe Mixing the noise with the car- 





rier wave results in a noise process with a PSD with half 





N 
the power density or Sm(f)=—P + which is the PSD at the in- 


put of the matched filter. The output of the matched filter 


has power spectral density 





N 
Sou (F) = Sn (FE (FY) =SNH(f) - (3.14) 
To integrate equation (3.14), Parseval's theorem [12] 


allows the time-domain impulse response function, h(t)=p,(t), 








to be used instead. Therefore, the power of the noise at 





the output of the matched filter is 


= (3:5) 


Since p,(t)=1 for O<t<TJ, and has zero value otherwise, 


equation (3.15) resolves to the noise power at matched fil- 


ter output. 
5. Reference Receiver Error Performance 


Using the output of the matched filter, the signal of 





interest and nois nter th decision block wher th 1g 








ceiver makes bit determinations. Using the signal of inter- 





A 
est input, W(KT, =F day the decision statistic is a Gaus- 
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Ad, 1, No, 


; ie This is used to de- 


Sian random variable, v( 





termine the probability that the decision block output bit, 


A 


da, LS in @rror when the sent bit is .d.,=]=—l.- Daviding by 





N.T, 
the standard deviation of the decision statistic, ira 





one finds the area under the right-hand tail of the Gaussian 


distribution, resulting in 


Pr(d,_,=11d,, =-1)=Pr(y(kZ,)>0) 
= AJT, (3.16) 


VNo 


Using the equivalence A7,=2E,, where E, is the energy 





per bit, the resulting probability of bit error is identical 








to BPSK, as expected since the signal is antipodal [7]. 


P (d =lld =-l)=Q 
r dec k-1 
No 


2E, 
-of a 


The theoretical error performance for an antipodal sys- 





(317) 


tem plotted using a range of signal to noise (SNR) ratios or 


E 
— (bit energy to noise power spectral density) ranging 
0 


from 0 dB to 10 dB is shown in Figure 3. This curve repre- 





sents the theoretical best error performance of the refer- 





nce receiver with uncoded BPSK modulation and was generated 


using MathCAD (v.14.0.0.163). 





; v 2 : : . : : 
1 The notation N(X,o ) is used in this thesis to describe a random 
variable as a normal, i.e., Gaussian, random variable with mean X and 
variance oO . 
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Bit Error Rate (BER) 


Figure 3. 


E/N, (dB) 





E, 
Range of — 


0 


B. MODELING THE REFERENCE RECEIVER 


Th 





Simulink. 


variabl 


By setting this variable to three, 
100,000 bits for each SNR level over th 


dB. The AWGN channel’s calculations wer 


Le 





each bit. 





Optimum Bit Error Rate for a BPSK System over a 


referenc receiver was simulated with MATLAB and 


The MATLAB file, shown in Appendix A, 


“sim_length” to control number of bits 


used 


simulat 
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the 





ce. 


the simulation ran for 
range of 0 dB to 9 


independent for 


‘i Simulation Results 


The reference receiver’s simulation results are shown 





in Figure 4. The black line represents the BPSK theoretical 
bit error rate performance, and is the same curve from 


Figure 3. The simulation bit error rate results are the red 





circles on the plot, representing the bit error rate (BER) 
for each SNR level. The calculated error rate was the total 


number of errors divided by the total number of bits sent. 





The results of the simulation matched the general trend of 


the reference curve, indicating an accurate model. 





Reference Receiver Performance in AVVGN 


35 BPSK Theoretical 
© Reference Receiver 


BER 








E/N, (dB) 





Figure 4. Reference Receiver Performance in AWGN Over a 
Range of SNR Values and Compared to BPSK 
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2. Transmitter Simulation Design 


The reference transmitter simulation generated a single 








user’s Bernoulli binary random data stream at the bit rate 
and a 15-stage shift register maximal length PN sequence at 


the chip rate. The simulation mapped unipolar binary out-— 











puts to antipodal binary before mixing. The resulting 
chipped signal mixed with the carrier wave, a cosine func-— 
tion block with a carrier frequency of 1.2288 MHz. The re- 


sulting signal then passed through a fixed gain before en- 





tering the AWGN channel block. The simulation’s parameters 





are listed in Table 1 and are based on IS-95 parameters [5]. 


Table 1. List of Reference Transmitter Simulation 























Parameters 
Parameter Set Value 
Chip Frequency R, 1.2288: MHZ 
Bit Frequency R, AO 2 se 
Carrier Frequency /f, 4.915200 MHz 
Sample Time T, 1.27E-08 seconds 
User Gain A 1 (no units) 














Simulink’s AWGN block provided the AWGN channel simula- 
tion using the parameters shown in Table 2, which the MATLAB 


file set prior to starting Simulink. 
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Table 2. List of AWGN Channel Simulation Parameters 

















Parameter Set Value 
Initial Seed 45 (default setting) 
Mode Signal to noise ratio 
E,/N, (dB) 0 to 9 dB (set by MATLAB 
file) 








Number of bits per symbol 





Symbol Period 1 
R, 














The original data goes to the workspace for use in the 


receiver and error calculations. The combined signal then 





mixes with a cosine wave function with a carrier frequency 


that is four times the chip rate and sampling frequency that 





is 64 times the chip frequency. The user gain is unitary 


and has no effect on the model with the parameters as set. 








Integer relations between the chip and bit frequencies wer 


arbitrary, but we require f,>R, and R,>>f,. 


3. Reference Receiver Model 





The receiver shown in Figure 5 simulates the referenc 





receiver in Figure 2. It consists of a downconverter and 





filter subsystem, a baseband processor, and an integration 
block that compares the demodulated signal to a threshold to 


determine individual bits. The error rate calculation block 








compares the threshold output against the original data gen- 





erated by the Bernoulli binary sequence block in the trans-— 





mitter. However, this visual output is for the user and a 
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separate sink saves each bit to the workspace in a single 


vector for later error calculations and plotting. 


oY 


Error Rate 
Calculation 


cope 


1 > Corvert 
Original Data In 


TX Error Rate 


Original Data M, 
- »Rx Calculation 


Data Type Conversion 
errors 
> Reference RCVR_BER 


2 Pir(t) Sig out}—®Signalin + Decs Vat—® <= jntegrationT hreshold > To Workspace 
signal received = 
" Baseband pare > Data_Recovered 
, 5 CTS T nstant 
nae: Processor To Workspaca 
el 

















Figure 5. SIMULINK Reference Receiver Model 


The downcoverter and filter subsystem mixes the signal 
with a cosine wave with the same parameters as in the trans-— 
mitter. Next, the signal passes through a discrete-time 
low-pass filter (LPF) using a unitary pulse that is 64 sam- 
ples in length to shape the output. This simulates the ef- 
fect of the matched filter and eliminates the double fre- 
quency. Finally, the receiver downsamples the output with 
downsampling factor of 64 before proceeding to the baseband 


processor, which simulates the reference receiver’s sampler. 





The baseband processor models the de-spreading, matched 





filter, and sampler from Figure 2, and mixes the received 


signal with a PN sequence that matches the transmitter. The 





resulting signal passes through a matched filter of 64 ones 
and downsamples 64 times. The processed signal output then 
goes to the decision block to determine bits. The integra- 
tion block performs the threshold comparison and decision 
functions producing a unipolar binary signal output, and 


sends the results to the workspace for comparison. The re- 





ceiver simulation matched the theoretical performance, so it 
provides a point of reference for performance comparison 
with the cyclic-PN code receiver in the next chapter. 
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IV. CYCLIC PSEUDO-NOISE CODE RECEIVER 


The cyclic PN code receiver solves the scalability 


problem while maintaining the BPSK error performance of the 














reference receiver. Its use of DFTs implemented through the 











FFT algorithm eliminates the need for multiple matched fil- 











ters to extract each transmitter’s data. Rather, a single 








receiver block performs all the same functions for all the 








transmitters except for bit decisions. 


A. CYCLIC RECEIVER DESIGN 


The cyclic PN-code design uses two spreading codes, 





with each transmitter using a cyclically rotated version of 
the same short code mixed with a longer common PN code. The 
shifted code repeats for each bit period and allows the re- 
ceiver to identify the transmitter and demodulate data. 


When mixed with the long PN code, whose duration is much 





greater than a single bit, the resulting pseudo-noise s 


quence is unique and improves the LPE properties of the sys- 





tem. The DFTs allow a single receiver to analyze all trans- 





mitted signals simultaneously and determines the number of 








shifts of the short PN code in each signal as well as the 








data content. From this, the receiver determines transmit-— 
ter identity and makes bit decisions. The theoretical bit 
error rate performance of this receiver is the same as the 





2E, 
reference receiver, P,=Q| |—*|. 
No 
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1. Cyclic PN-Code Transmitter 


The cyclically shifted PN code transmitter mixes the 


long PN code, bit) ¢ with a cyclically shifted shorter code, 





c(t.7(p)) , where t is time, vy is the length of the shift, 
and p is the transmitter number. The transmitter number 


and quantity of shift are related by an integer factor k, 


resulting in a relationship of y=pk. For this analysis, 





k=l, so the shift number will be the same as the transmit- 
ter’s number, y=p. The shifts in c(t,7(p)) are not calcu- 


lated by each transmitter but are preprogrammed prior to 
system emplacement. The block diagram for the transmitter 


is shown in Figure 6. 





th 


Figure 6. p Cyclic PN Code Transmitter 





The cyclically shifted PN code repeats within each bit, 


allowing for demodulation and identification. It can be 








N-| 


represented as e(t.y(p))= cl i-7(P)], Pr (t-it.), where p is the 


i=0 


transmitter number, t is time, y(p) is the number of shifts 





of the reference PN code for transmitter p, and p,(t)=1 if 





O<t<T. and Pr, (t}=0 otherwise. The sequence cli 1\y 2. OO= 


24 


tained by cyclically shifting the N chip sequence, cli], y 





shifts to the -right... Hach cli] ¢{-1,1} and p does not have to 


equal vy. 


The longer PN code, b(t), provides security by not re- 


peating for the duration of the message. The two PN codes 





are generated and mixed at the chip rate R.,. Using two 


codes creates a unique spreading sequence for each transmit-— 


ter. 


The resulting unique PN code then mixes with the an- 


tipodal data signal m, before upconversion to a BPSK trans-— 


mission at a predetermined carrier frequency f.. The an- 








tipodal data signal can be represented as m,(t)=>d,p, (ti), 


where d, is the i” bit and p,,(t-iT,) is the pulse function 


fon. thes f° ble ames Finally, the transmitter amplifies the 





Signal by a pre-determined factor A before transmission, 


yielding the transmitted signal 
s(t)=Ab(t)c(t,y)cos(27f,t) . (4.1) 


2. Cyclic PN Code Receiver 


The signal in equation (4.1) travels through an addi- 





tive white Gaussian noise (AWGN) channel and arrives at the 








receiver depicted in Figure 7 with the form shown in equa- 


tion (4.2). 


r(1) = AL de(1.7)(s)eos(2a,.2) +m(0) (4.2) 


The signal enters the downcoverter block and changes to 
the form shown in equation (4.3). 


2D) 


= AY d,p, (t-iT, \c,(t,7)b(t)cos?(2nft)+n(t)cos(2rf.t) (4.3) 


s(t)+n/(f) b(t) C(k,0) 


DFT Le) —+[ ETL 


cos( 27 ft) 





Figure 7. Cyclic PN-Code Receiver 


Based on the trigonometric identity, cos" (x) =5(1cos(2x)) 


[ll], this changes to the form shown in equation (4.4), with 


the signal received now in three parts the desired signal 





term, the double frequency term, and the noise term. 


roe 5 aire (1 iT, )c(t, 7) b( +S Darl iT, )c(t, 7) b(t) cos(4zf.t) 


+n(t)cos(2z ft) 


(4.4) 


Next, examining the desired signal term allows deriva- 
tion of the mean of the decision statistic. Noise will be 


examined separately in the next subsection. Mixing with the 





long PN cod ffectively cancels the factor b(t) from the 


first two terms on the right hand side of equation (4.4). 





The LPF then removes the double frequency term, giving the 


output 


=F Da LevAlle-( (j+1)T,)/T, | (4.5) 


where A(x) =1-|x] if |x| <1 and A(x)=0 otherwise. 
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The receiver then samples the signal f (t) at the chip 








rate prior to it entering the DFT in order to analyze a sin- 











gle bit. 
A 

g[n|=>deln—7], (4.6) 
Analysis of the signal-only term leads to the mean of 
the decision statistic. Starting from the output of the in- 
verse DFT (IDFT), the output of the receiver is shown in eq- 
uation (4.7). C(k,0) represents the complex conjugate of 
the. “DET Of 9 the. short “PN code “with, Ta -shufie, Prem 





calculating this result and storing it in memory rather than 





performing a real-time DFT and complex multiplication helps 


improve scalability. 





i[n] = 1DFT| {C(k,0) | DFT(s[n]) | (4.7) 





Multiplying the two DFT outputs is equivalent to circu- 





lar convolution in the time domain [12]. The operator ® 


indicates a circular convolution. 
H (k)X (k) = DFT {h[n]® x[n]} (4.8) 


The complex conjugate, indicated by *, in the fre- 





quency domain is equivalent to the time-reversed signal in 








the time domain [12]. In this instance, the complex conju- 


gate for the discrete time domain signal in reverse is real, 








so the complex notation drops out. 
X\k 
aS [| (4.9) 
x[-n]<o X°[k] 


A circular convolution is the same as an ordinary dis- 








crete time convolution, except that the time-shifted signal 


is circularly shifted, vice an ordinary shift [12]. 
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n{n]® x(n] = h[m]x[n—m], (4.10) 


m=0 


Applying equations (4.8), (4.9), and (4.10) to equation 








Ad,NO, , 
(4.7) yields an IDFT output of ie as shown in 
equation (4.11). This: holds for 0<m=N-l1, where o,, is the 





Kronecker delta and it is assumed that c[n] has low correla- 


tion with c[n-i],, for all O<i<N-1. The square of a shifted 


PN codes is one and the sum of that result over N-terms is 





No. 
i[n] = IDFT| {C(k,0) | DET (g [*}) | 
= g[n]®c[-n] 
= Seln—7), Sel] a 
-hYelm-r]yelm—], 
~ AGNS 7 
y) 
When n#¥y, using a circular convolution definition 
with a Kronecker delta gives a result of l(n#¥y)=0.  Choos- 


N-1 
ing Yic[m-y], c[m—n], ~ N6,,, allows the result shown in equa- 


m=0 





tion (4.11). In the case of multiple transmitters’ signals 


arriving simultaneously, I[n] indicates whether transmitt 





nm is transmitting (if I [n]| is large) and the value of 





Ler 





the 


sent data bit (sgn(J[n])) . Thus I[n] for O<n<N-1 simultane- 





ously indicates which of the N transmitters are active and 


the bits they have transmitted. 
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3. Noise Calculations 


Noise at r(t) is additive, white and Gaussian and has 


Ny 
2 





u=0, oOo =+0, and a PSD of The downconversion 





N 
changes the PSD to 1 due to the average power of the 
squared cosine function, which is ie 


After the LPF, the variance of the noise is no longer 


infinite, but limited by the bandwidth of the filter. 


of No : 
a=] SH) 


Ny +00 
> a 
_ NRT, 
4 
_NjR,N 
oT a 


h(t)| dt 





(4.12) 


where H(f) is the frequency response of the LPF and sam- 





pling does not affect white Gaussian noise. Therefore, th 





noise at g(r], the input to the DFT, is still a Gaussian 


NOR 
4 


random variable with Zero mean and or 








N,R,N 
y[o Ae), The samples, gln], are independent because th 


input to the LPF is white noise, the LPF’s impulse response 


has a duration equal to the chip duration, and the LPF’s 





sample rate is equal to the chip rate. Therefore: 
E{g[m]s[p]} = 


NR,NO,, p 


4.13 
i ( ) 
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The output of 


the mixer in the frequency domain is: 


L(k)=G(k)C(kI0) . (4.14) 








The expected output of the receiver is then the IDFT of 


L(k). The product 





main is equivalent 


main [12]. 


The expected 


this function is: 


The variance 


of the two functions in the frequency do- 


to a circular convolution in the time do- 


I[n]=IDFT {G(k)C(k10) | 


= g|[n]®@c[-n| (4.15) 


value or mean of the noise component of 


N 


E{I[n}} 21S elmeln-n], | 


N-l 


= E{g[m]}c[m-n], (4.16) 


of this function is the variance of the 


receiver output and therefore the variance of the signal 





(which allows probability of error calculations), 


o£) S¢(m)elm—n), 2e(r}e(n~n)} 


m=0 p=0 


E{g(m) \[e(m=n), | (4.17) 
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4. Theoretical Error Performance 


At the output of the IDFT block, the signal’s decision 








desde sf Ad,N N,R,N* ‘ae 
statistic is I[n]~N 5 ; 5 . This is used to find the 








2E, 


probability of bit error, of N 
0 


} the same performance as 


a BPSK system. Using the mean and variance calculated above 





in equations (4.11) and (4.17), and the equivalence 


E,=0.5A°T,, the bit error rate is 


P, =Pr(I(n)>0ld, =-1) 








AN AN 
jee 
a (4.18) 
ae pan 
4 4 


2E 
No 
B. MODELING THE CYCLIC PN-CODE RECEIVER 


The proposed receiver was simulated in MATLAB and 





tested for 100,000 bits at each of a range of signal to 





noise ratios Ga from 0 dB to 9 dB. A single transmit-— 
0 


ter was simulated with a one chip cyclic shift (y=1) in the 





base short PN-code. The simulation used an AWGN channel and 
BPSK modulation for the carrier. Perfect phase synchroniza- 


tion was assumed. The MATLAB code for the simulation is in 





Appendix B. 
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ce Simulation Results 


The cyclic PN code receiver simulation results are 
shown in Figure 8. The black line represents the BPSK theo- 


retical bit error rate, and is the same curve from Figure 3. 








The simulation bit error rate results are the blue squares 
on the plot, representing the BER for each SNR level. As 


with the reference receiver, the calculated error rate was 





the total number of errors divided by the total number of 





bits sent. The results of the simulation matched the gen 
eral trend of the analytically derived curve, indicating 


that model’s accuracy. 








Figure 9 is an example of the IDFT output for a single 
bit from a single transmitter with a single shift of the 
reference short PN code, i.e., y=l. The signal from the 


transmitter of interest stands out clearly from the 64 pos- 





sible signals received during this bit time. 


Cyclic PN-Code Receiver Performance in AVWWGN, Eb=0.5 
= =e tT 


es z 

















ssocces t--------J--------------------------- 
BPSK Theoretical 
O Cyclic-PN Code Receiver 











BER 




















i) 1 2 3 4 5 6 it 8 9 10 
EWN, (dB) 


Figure 8. Cyclic PN Code Receiver Performance in AWGN over a 
Range of SNR 
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IDFT output for bit #1, SNR=9 dB 
————SS ———— ee = ee aap 


6000 + 4 


4000 + 


2000 ~ 


IDFT output value 
Oo 
i 


-2000 + 4 


4000+ 4 





-6000 + 4 











0 10 20 30 40 50 60 
Discrete Time n, y=1 








Figure 9. IDFT output for a single bit duration, at 


Ef = 
if, =9 dB with y=l 


Figure 8 shows good concurrence between simulation and 





theory. Figure 9 shows that I[n] at the receiver indicates 


reception of a +1 from transmitter number one and only noise 


from the other transmitters. A mechanism is still needed to 





distinguish an active transmitter’s signal from noise only. 
This is not addressed in this thesis except to note that 
conventional methods such as uSing start and stop bit pat- 
terns or employing two thresholds to make “+1”, “not trans-— 


mitted” and “-1” can be used to address this issue [13]. 
2. Transmitter Simulation Design 


The simulation of the cyclic PN code transmitter used a 





fixed R,=l and then generated other variables from its mul- 
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tiples. Unlik th referenc receiver, these differences 


were relative and not on an absolute scale. The simulated 





transmitter mixed the fixed data sequence with two separate 


PN codes, one of which was the cyclic rotation of the short 








PN code. The resulting signal was sampled at the chip rate 





R,.=64R, and then upsampled with upsampling factor of 256 be- 
fore mixing with the cosine vector at frequency, 
f. =16R, =1024R,, and a fixed amplitude. The transmitter pa- 


rameters are in 0. 


Table 3. Parameters for Cyclic PN Code Transmitter 
































Simulation 
Parameter Value 
KR, Abit rate) 
R. (chip rate) 64R, 
f, (carrier frequency) 16R, =1024R, 
R, (modulated signal sample rate) 16f. 
T. (sample time) VR 
Data bits 50,000 





The cyclic PN code simulation did not use the AWGN 











block in SIMULINK. Rather it generated noise as a vector of 





the same length as a single upsampled bit based on the SNR 








LOT a particular loop and the noise PSD using 
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R.N 








Noise=  Tandn(1,R,). The noise PSD was then based on the 
; E, 
energy per bit by No =~ sua * with if, fixed. at: 0.5% The 
10 10 


MATLAB function “randn” generated a vector of random numbers 


with length equal to the sampling rate with a coefficient 





based on the noise PSD, allowing variable noise based on SNR 
and vector addition. The effect on the transmitted signal 


was equal to the effect of AWGN. 


3. Cyclic PN Code Receiver Model 


The cyclic PN code simulation downconverted the re- 
ceived signal with noise by multiplying the data vector by 
the same cosine vector used in the transmitter. To de- 
spread the signal with the long PN code, the simulation up- 
sampled the PN code before mixing. A square pulse whose 


length was the chip duration performed the low pass filter 





function through convolution with the de-mixed signal. The 
simulation sampled the LPF output at the chipping rate, by 
downsampling using the ratio between sampling and chipping 


rates. 


MATLAB’s Fastest Fourier Transform in the West (FFTW) 





algorithm, the “fft” function, acted as the DFT block pro- 
posed in the theoretical model. The simulation directly 


calculated the DFTs of the signal received and reference 





short PN code, the complex conjugate of the short PN code, 





and mixed the DFT outputs. The MATLAB function “ifft” acted 








as the IDFT block, and a loop made bit decisions based on 





the IDFT output. 
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The cyclic-PN code receiver simulation performed well 





and matched the theoretical BPSK error rate. In addition to 





the cyclic-PN code receiver, a coding mask receiver was pro- 





posed as a possible alternative. The next chapter analyzes 





the coding mask receiver and compares its theoretical error 


performance to the reference receiver. 
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Vv. CODING MASK RECEIVER 





As an additional point of comparison, one possible al- 





ternative to the cyclic PN code receiver is a coding mask 





receiver. This design relies on similarities in the trans- 








mitters’ PN codes to perform a two-stage detection and de- 


modulation. This receiver is not a viable solution as a 





system if many transmitters are expected to be simultane- 








ously active as it does not improve scalability in this case 





and it has a greater error rate than the reference receiver. 








However, it requires fewer processes for single transmitters 





and could be a viable alternative design. 





A. RECEIVER DESIGN 


The coding mask receiver subdivides the full body of 


transmitters into groups, which will be called subsets for 








the remainder of this thesis. Each subset has related PN 


codes with some common chips to speed-up processing and this 








is the “mask.” The receiver’s signal detection circuit mix- 





es the combined signal received with the base PN code of its 





subset and then its mask to determine signal presence. When 





a Signal from a particular subset is detected, all receivers 








within that subset process in parallel to extract the signal 
received. This design has non-linear growth and while it is 
of similar complexity, efficiency is gained through savings 
in total number of processes. This savings would be most 
beneficial in a software defined radio implementation since 
the savings in processing is realized by simpler and faster 


algorithms and reduced requirements on the microprocessor, 





Sid 





digital signal processor, or field programmable gate array 


(FPGA) responsible for the signal processing [14]. 


The coding mask receiver system would use the same 


transmitter as the reference receiver as shown in Figure l 





and generate a similar transmitted signal in the form 


s(t)=Am, (t)c, (t)cos(27f.t) . Here, m, (t) is the antipodal data 





with the subscript p indicating the transmitter number and 
having the same definition as equation (3.1). The PN code, 
ealtyy is unique and antipodal for each transmitter and the 
subscript again indicates the transmitter number. The local 


oscillator, cos(2zf.t), is also BPSK, and the coding mask 





transmitter also amplifies the transmitted signal by a pre- 
set gain A. 

The signal then passes through an additive white Gaus-— 
Sian noise (AWGN) channel, which adds the noise term n(t) to 
the combined signal. The signal received is, therefore, 
r(t)=Am, (t)c, (t)cos(2af,t)+n(t) . This signal ‘then -enters the 


coding mask receiver shown in Figure 10. 
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Figure 10. Code Masking Receiver Diagram 





In Figure 10, a receiver for a single subset is shown. 





The PN codes of this subset have some common points allowing 


the mask to be pre-calculated and preprogrammed into the de- 


tection branch. The first PN code in this subset is Colt) s 


where the subscript O indicates the ordinal numbering of 





the transmitters in the subset. The first stage accom- 


plishes signal detection by mixing the downconverted signal 
with the first PN code in the subset mea the carrier 
cos(27f.t), and the mask k(t) before entering the matched fil- 


ter. The mask represents the common points of all PN codes 





in a subset, with 1 used to indicate they are the same and 
0 used if even a signal code is different at a particular 


chip. The results is k(t) € {1,0} and not e {1-1} . The signal 





received at this point is similar as in the reference re 





ceiver with three components, a signal of interest, a double 


frequency, and noise [7]. 
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r'(t)=—m(t)k(t)——m(t)k(t)cos(4zf.t)+n(t)c,(t)k(t)cos(27f.t) (5.1) 
The matched filter’s impulse response h(t) is a pulse 


function p,,(t-iT,), where the value of the function is 1 


during the bit interval and 0 otherwise. The resulting out- 








put is the signal component of interest output [7]. 
m(t)k(t)*h(t) 

© (5.2) 
=5 J m(a) )h(t-a)da 


However, the integrand is only non-zero when t-T,<a<t, 





changing the integration interval. A substitution can be 





made for the modulated signal, whereby the antipodal data, 





m(t), can be represented as a single bit at a discrete time 
instance rather than as the sum of bits and the pulse func-— 


tion, as was shown earlier: J= id pat a-iT,)=d,,. Addi- 





tionally, the impulse response of the matched filter taken 


at a discrete time instance where t=f7T, results in the sub- 





SeteuL LON A(IT, -a) = p, (IT, -@) . Using these two substitutions 


changes, the matched filter output in the detection branch 


to a different form of the mean of the signal received at 








that point. 
A t 
y()=> J m(a)k(a)h(t-a)da 

iF, 

A 

y(IT,)=— | d,,k(a) p,(41,-a)da (353) 
2 (yn, 
A 
= ay Ayaft T, 
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This result is identical to th referenc receiver’s 





signal of interest with the exception of f#, the fraction of 
T, for which the mask is one, i.e., B= — k(t)dt . Since the 


matched filter is a low pass filter, the double frequency 


term is filtered out. 

















N 
The PSD of the noise at the receive antenna is ke 
As in the reference receiver, when this noise is mixed with 
: : 1 ' Ny 
the carrier, the resulting PSD is ain he However, the mask 
has two possible states, k(t) =0 or k(t)=1. In the former 








case, the resulting mean (f#), variance CE Ss and PSD (S$) 





of the noise would be zero. In the latter case, w=0, 


2 N ; ; 
Oo =+0, and ae which is white Gaussian noise. However, 


to account for the fraction of the mask that equals 1 ina 





given period of 7J,, the term # can likewise be used as in 





the signal of interest, giving a resulting input PSD of 


ea 
4 


in 





The noise PSD at the output of the matched filter 





2 NP 


is therefore Sou (£)=|F (F) 4 


Its integral is the noise 


power at the matched filter output 





Oe j la (f)) ae df - 7 f JH (Ff) af . Applying Parseval’s theo- 


N 
rem [12], changes this to oF =O | he (t)ar . And this in turn 


solves to the noise power or variance 
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9 NBT 


5.4 
4 ( ) 


B. PERFORMANCE ANALYSIS 


Using the sampled matched filter output for the signal 
of interest and noise, the decision statistic for signal de- 


Ad,_,BT, NPT, 
Dt 








tection is v( ). The probability of bit error 


is the same as the probability of deciding a “1” was sent 


when a “-1” was sent. As shown in equation (5.5), the prob- 





ability of error for a “-1” sent is equivalent to the prob- 








ability that the mean of the signal plus the noise N is 


greater than zero. 


Pr(d,,, =11d,, =—1) =r AF w >0) 


ar (3:35:) 
lf 


In terms of signal power required, the best possible 








condition would be f=1, when all PN codes ar th sam 





However, this condition would make multiple access impossi- 





ble-the receiver would be unable to distinguish between the 





signals. 
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Error Performance for SNR = 5 dB with variable Beta 





Oo 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Beta 


Figure 11. BPSK Error Performance for Fixed SNR 
and variable PB 
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Effects of Beta on BER 


BER 


| BPSK Theoretical |.:. 
} —&— Small Beta, 5/64 
—— Half Beta, 0.5 
||? Large Beta, 58/64 |: ; 
i) 1 2 3 4 5 6 rd B g 10 
E,/N, (dB) 





Figure 12. Comparison of Theoretical BPSK and 
Three Values of PB 


Fewer chips in common only increases the bit energy re- 





quired to achieve a similar bit error rate as the reference 


receiver. The resulting system probability of bit error is 





worse than reference receiver for the same signal to noise 
ratio. The first Q function is the error of the detection 
branch where # represents the effect of the mask. The sec- 


ond Q function is the demodulation error of the receiver 


branches in the subset. The result is a performance worse 
2E 

than 20) /— |. 
Ny 


ros) =o 288 ).of 25. (5.6) 


This receiver was not simulated due to time constraints 


and its significantly poorer BER performance as compared to 











the receiver in Chapter IV. Additionally because of its in- 





ferior theoretical error performance, the overall system 








performance was not modeled and compared to the reference or 





cyclic PN code receivers. Chapter VI, therefore, does not 








include the mask receiver’s performance in its comparison 


and analysis. 
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VI. RESULTS AND COMPARISONS 


The reference receiver and cyclic PN code simulations 





ran were of the same duration, 100,000 bits per dB level. 





Their results both closely matched the theoretical un-coded 
BPSK performance in an AWGN channel, as expected. Of note, 
both simulations did not address fading or synchronization 


issues. 


A. BIT ERROR RATE PERFORMANCE 





The reference receiver and cyclic PN code DFT receiver 





have identical BER performance, as shown in the analysis in 























Chapters and IV, and confirmed in the simulation results 





shown in Figure 13. In both cases, a SNR of approximately 7 





dB was required to achieve an error rate of 1 in 1000 and 





approximately 8.5 dB for 1 in 10000. 
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Comparing Receiver Performance in AVVGN 


BPSK Theoretical 
© Reference Receiver 
O Cyclic-PN Code Receiver 


BER 





o 1 2 3 4 65 6 7 8 93 10 
E,/N, (dB) 





Figure 13. Comparison of BPSK Theoretical (curve), Reference 
Receiver (circles), and Cyclic PN Code Receiver (squares) 


B. SCALABILITY COMPARISON 


Examining the reference receiver system with N trans- 








mitters, it initially processes the received signals with a 





common antenna, low-noise amplifier, local oscillator, and 








low pass filter. From there, such a receiver requires a 
unique matched filter, sampler, and bit decision block for 


each transmitter in the system. 





The cyclic PN code design in Figure 7 and the reference 


receiver in Figure 2 have similar components from antenna to 








LPF. Looking at the number of multiplications required al 
lows for a scalability comparison. For the reference re 
ceiver, there are N branches from the LPF forward. Each 


branch must analyze N signals, with despreading, filtering, 
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sampling, and bit decisions. Therefore, a good estimation 


of the scalability is total number of multiplications re- 











quired for despreading in all the branches. Sign multipli- 
cations (recall c[n]e{-L1}) are not as complicated as other 


multiplications, so a coefficient, a<l, is used to indicate 


the ratio of the complexity of a sign multiplication to the 





complexity of a complex number multiplication. For each 





branch in the linear-growth receiver, there are N sign mul- 





tiplications per bit duration because there are N chips per 








bit. Therefore, the total effective number of complex mul- 





tiplies per bit duration is 
C =N(aN) 
=N’a 


This value is accepted as the measure of complexity of the 


reference 


(6.1) 





reference receiver. 





For the cyclic PN code receiver operating in the same 





scenario, it must complete N sign multiplications per bit 


duration to despread the long PN code, N complex multiplies 











in the last mixer, 1 DFT, and 1 IDFT. The FFT and IFFT each 








require O5Nlog,N complex multiplies [6]. Combining these 


gives the total number of complex multiplies per bit dura- 





tion for the proposed receiver, which is accepted as its 


measure of complexity. 


Coctic-py = AN + Nflog, N+N+ Nflog, N 








(6.2) 
= N (log, N+1+@) 
The ratio of the complexity of the reference receiver 
to the complexity of the cyclic PN code receiver is 
Na 





, where the value of @ depends on the bit resolu- 
log, N+l+a@ 
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tion of the numbers multiplied. The complexity of one com-— 
plex multiply with real and imaginary parts each represented 


with r bits is equivalent to the complexity of 4r sign mul- 





tiplications. If the number of transmitters exceeds ap- 


proximately 300, the proposed receiver is a more efficient 
implementation, assuming a=\, which corresponds to eight 


bit resolution. The relative growth in the size of the two 


systems is shown in Figure 14. 


x 10° Comparing Receiver Scalability 


| Reference Receiver | 
: | — — - Cyclic-PN Code Receiver |! 


Number of Complex Multiplication Equivalents 





0 100 200 300 §=©6©400 500 600 700 800 goo =. 1000 
Number of Transmitters 


Figure 14. Scalability Comparison for N#=1000 and a=, 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


As can be seen in Figure 8 and Figure 13, the proposed 


design matches theoretical uncoded BPSK performance (which 








also matches the reference receiver). The IDFT output, as 











shown in Figure 9 for a single transmitter, clearly indi- 





cates the shift number and bit value. The simulation used 


an N=64-chip PN code, allowing up to 64 transmitters in 





the model system. For a full system, 64-bit decisions would 





be required for each 7, interval. N=64 was chosen for il- 
lustration purposes only, and much larger values are practi- 
cal. The cyclic-PN code receiver maintains the BER while 


improving scalability, making it a viable solution. 





A. PERFORMANCE COMPARISON 


1. Reference vs. Cyclic PN Code 


Since the receivers had similar error performance for 





an uncoded signal in AWGN, it stands that the key difference 





is receiver complexity. In this, the cyclic PN code re- 


ceiver has superior scalability when the number of transmit- 





ters is large, as shown in equations (6.1) and (6.2), as 





well as in Figure 14. 


2. Reference vs. Code Mask 


The theoretical bit error rate performance of the code 





mask receiver was inferior to th referenc receiver. The 








code mask receiver’s advantage lies in the fewer required 





processes to receive a single signal, but its higher error 





rate offsets this advantage. 
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3°. Code Mask vs. Cyclic PN Code 


The cyclic PN code receiver has the same bit error rate 





performance as th referenc receiver, but superior scal- 


ability and, therefore, is the better choice when the number 








of transmitters is large. Additionally, the DFT-based fil- 








tering in the cyclic PN code receiver was more efficient in 





terms of space, but potentially less efficient for a single 








signal as the total number of processes required is the same 


regardless of the number of transmitters in the system. 





However, for a large number of transmitters, the coding mask 


receiver is inefficient, as it requires two stages of proc- 





esses, resulting in more overall complexity than either the 
reference or cyclic PN code receivers. The coding mask re- 
ceiver’s theoretical BER, even under the best conditions 


would be twice that of the cyclic PN code receiver. 


B. FUTURE WORK 


1. Improving the Cyclic PN Code Model 





The models presented were simulated in an AWGN environ- 





ment. Simulating the same models in a fading and jamming 


environment would better match potential real-world scenar- 








ios. Large networks of sensors deployed in a mountainous or 
urban environment would suffer fading or jamming effects in 
addition to AWGN. Validation of the cyclic PN code re- 
ceiver’s error performance in these environments is impor- 


tant. 


MATLAB’s “randint” function and Simulink’s equivalent 





block generated the simulations’ random unipolar binary se- 











quences. Future researchers could develop a truly orthogo- 
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nal cyclically shifted PN code and validate its performance, 





especially the prevention of inter-transmitter interference. 





Additionally, future research should synthesize an effectiv 








forward error correction code using BPSK or QPSK for the cy- 


clic PN code receiver. 


Quantifying the LPE features expands the validity of 
the proposed communications system. This would require a 
larger model with multiple transmitters and a non- 


cooperative interceptor. Additionally, models of the most 





likely operating environments, which account for potential 














adversaries, are necessary in order to analyze potential 


system performance. 


This thesis’ models also need a synchronization analy- 








sis. Both the reference and cyclic PN receiver in this the- 
sis relied on assumed synchronization to make bit decisions. 
Synchronization blocks or a non-coherent form of modulation 


might solve this problem. 


2. Improving LPI and Data Rate Performance 


Strong encryption would improve the systems defenses 


against exploitation of the intercepted signal. However, 





encryption would not improve overall LPE as adversaries 


could potentially exploit signal presence and geolocation 








through direction finding to localize transmitters. 


The systems designed and tested in this thesis used 
BPSK modulation. Modulation techniques like QPSK, M-PSK, M- 
frequency shift keying (MFSK), and M-quadrature amplitude 
modulation (MQAM) were not tested. Other modulation tech- 


niques could achieve higher data rates with trade offs in 





transmitter power, error rate, and exploitation. 
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APPENDIX A. REFERENCE RECEIVER SIMULATION CODE 


All code in this appendix was generated using MATLAB 





v7.6 (R2008a). The SIMULINK file that supports this simula- 











tion is provided via electronic media. 


ole 


Name: LCDR Frank Cowan 

Title: Reference Transmitter and Receiver 
Simulation Project: NPS EE Thesis 

Date Created: 22 Jan 2009 

MATLAB/Simulink Version: 7.6.0, R2008a 


AP ol? ol? 


ol? 


\o 


+ This file provides the variable inputs for a simulation of the 
reference transmitter and receiver for my thesis. All data is saved 
to files and will be plotted using a separate m-fil 





GLE $ clears the command window 


[3 


clear all; % clears workspace 








fe) 


% Transmitter parameters 

Chip_Freq = 1.2288e6; % Rc, set at 1.2288 MHz 

Bit_Freq = Chip_Freq/64; % Rb, set at 19200 Hz, 1/64 of Re 
Carrier_Freq = 4*Chip_Freq; % 4915200 Hz, frequency of carrier 





CarrierSampleTime = 1/(16*Carrier_Freq); %16 samples/cycle 
TransmitterGain = 1; % A, for Arbitrary gain 


SAWGN channel parameters 

EbNo_array = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]; % the EbNo increments 
InputSignalPower = (TransmitterGain*’2)/2; % used in the AWGN channel 
block for input signal power 

AWGN_bps = 1; % used in the AWGN channel block for bps setting 


AWGN_symbol_period = 1/Bit_Freg; % sets the symbol period for the 
AWGN in the channel 








% Reference Receiver parameters 
IntegrationThreshold = 0; % required parameter for the reference r 
ceiver decision block 











fe) 


% Conduct simulation 
Run_No = 1; % sets the Run_No to default setting of 1 


sim_length = 3; % same simulation length used for each loop - based 
on the time interval used, this works out to around 50,000 samples 
per loop 





for EbNo_idx = 1:10 

% Simulation 

Run_No % displays the run number in the work space to allow user 
to observe progress 

EbNo = EbDNo_array (EbNo_idx); % assigns the value of EbNo used in 


the AWGN block in the channel 
ao) 


sim('Reference_RCVR_Model_V1', sim_length); % calls the simula- 
tion and assigns the sim length 


° 


%6 Calculate the number of errors 
Bits_sent = length (Data_Transmitted)-3; % number of bits sent by 





transmitter. -3 is needed to account for delays in the system. 
Errors = sum(xor(Data_Transmitted, Data_Recovered))-1; % calcu- 


lates the number of bit errors by XOR'ing the data transmitted and 
data received. 

Error_Rate (EbNo_idx)=(Errors/Bits_sent); % calculates the error 
rate by dividing the number of errors by the total number of bits 
sent. This is the data that will be plotted after the loop. 








{23 


% save the number of bits created for the simulation run 

Simulation_Length (EbNo_idx) = Bits_sent; % saves the value of 
Bits-sent to an array for later analysis 

% increment counter 

EbNo_idx = EbNo_idx + 1; % increments the EbNo index and walks 
through the EbNo array 

Run_No = Run_No + 1; % increments the counter for the loop 
end 


save ('Error_Rate', 'Error_Rate') % Saves the Error Rate variable in- 
to a separate file 
save ('EbNo_array', 'EbNo_array') % Saves the EbNo_array variable in- 
to a separate file 








6 End of Reference RCVR 
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APPENDIX B. CYCLIC PN-CODE RECEIVER SIMULATION 


CODE 


All code in this appendix was generated using MATLAB 


v7.6 (R2008a). 


A. 


CYCLIC PN CODE RECEIVER SIMULATION MATLAB CODE 


Name: LCDR Frank Cowan 

Title: Cyclic PN receiver Code 

File name: PN_Code_RCVR_VF.m 
Simulation Project: NPS EE Thesis 

Date Created: 06 May 09 
MATLAB/Simulink Version: 7.6.0, R2008a 


AP oP lO AP AP 


ol? 


ol? 


Linked functions: upsample2.m 


% This files simulates a single transmitter system for the cyclic PN 
code 
% transmitter and receiver system detailed in my thesis. 


ele: % clears the command window 
clear all; % clears workspace 





% create data array 

Rb = 1; % bit rate 

Rc = 64*Rb; %& chip rate 

alpha = 16; % ratio of fc to Re 

fc = alpha*Rc; % carrier frequency = alpha factor times Rc 
beta = 16; % ratio of Rs to fe 


Rs = beta*fc; % sample rate 








Ts = 1/Rs; % sampling time - inverse of sampling rate 

N = 49999; % length of vector factor for non random data vector 

n = O:N; % length of non-random data vector 

Data_Array = (-1).*n; % alternating 1,-1 square wave —- not random 





Message_Length = length(Data_Array); % measures the length of the 
random data vector —- not used with the non-random data vector 


fe) 


Chipping_Matrix = ones(1,64); % 64 chips per bit - vector of 64 ones 








% create the unique PN code for transmitter on base code with one 
SHEE 

Short_PN = randint (1,length(Chipping_Matrix))*2-1; % base short PN 
code with no shifts 

Long_PN = randint (1, length (Chipping_Matrix))*2-1; % base long PN 
code - this really should be longer than one bit duration 

PN_Codel = circshift (Short_PN, [0,1]).*Long_PN; % creating the unique 
PN code for Transmitter #1s 


%& Energy per bit and amplitude calculations 
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Eb = 0.5; % energy per bit - this will remain unchanged each loop - 
final problem should use 0.5, but 2 seems to be working better 
A = sqrt (2*Eb*Rb); % peak amplitude 











% this is the variable noise loop 
SNR = [0,1,2,3,4,5,6,7,8,9]; % SNR in DB 





for h = l:length(SNR) % control for variable SNR loop 
No(h) = Eb/(10*(SNR(h)/10)); % No calculation - uncomment for 
noise loop 





& this is the transmitter/receiver simulation loop 
for k = l:length(Data_Array) % control for the simulation loop 
Chipped_Data = Data_Array(k) *Chipping_Matrix; % repeats the 
data bit Rc times 
Mixed_Data = Chipped_Data.*PN_Codel; % mixes the chipped da- 
ta with the combined PN code over one bit interval 
Sampled_Data = upsample2 (Mixed_Data,alpha*beta); % upsamples 


the data to the sample rate 





% build local oscillator 





time = (0: (length(Sampled_Data)-1))*Ts; % time vector 
Cosine_Array = cos (2*pi*fc*time); % cosine vector 


% build the modulated waveform that will be transmitted 
Modulated_Signal = A*Sampled_Data.*Cosine_Array; % mixes the 
cosine array with the sampled data - result is the modulated signal 


= prof Kragh changes. 


% simulate the AWGN channel and add noise to signal 
Noise = sqrt (Rs*No(h)/2)*randn(1,Rs); % Noise calculation - 
must be repeated for each bit in order to be independently random 


Received_Signal = Modulated_Signal + Noise; % noise is added 
to the signal 


% receive and process the signal 

Downcoverted_Signal = Received_Signal.*Cosine_Array; % mix 
with same cosine 

LongPN_up = upsample2 (Long_PN, (alpha*beta)); % upconverts 
the long PN code to mix with the data received at the sample rate 

Demixed_Data = LongPN_up.*Downcoverted_Signal; % mixes the 
data received prior to LPF input with the upsampled long PN code, 


Diey 











p = ones(1,Rs/Rc)*Rc/Rs; % averaging factor for downsampling 
from Rs to Rc - change fc back to Rc later 

LPF_out = conv (Demixed_Data,p); % evenly spaces convolution 
and averaging - final LPF output, must be downsampled still 

DFT_input = 
LPF_out ( (0+ (alpha*beta)): (alpha*beta) :64* (alpha*beta)); % Downsam-— 
pling of the LPF output from sample rate to chip rate 

G = fft(DFT_input); % DFT of signal received 
C = fft(Short_PN); % Reference PN code DFT 
C 
L 





0 = conj(C); % conjugate of reference PN code 
= CO.*G; % mix DFT outputs 
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l(k,:) = ifft(L); % IDFT of mixer output - matrix of data 
end % ends the transmitter/receiver system simulation loop 


$---back in the outer loop 


° 


%& make bit decisions from IDFT output - idft_out_x_axis 
for m=1:length (Data_Array) 
if 1l(m,2) >= 0 


Bit_Decision(m) = 1; 
elseif l1(m,2) < 0 
Bit_Decision(m) = -1; 


end 
end % ends the bit decision loop 





° 


% convert both in and out data to unipolar binary 
Data_In = (Data_Array + 1)/2; 
Data_Out = (Bit_Decision + 1)/2; 
% determine error rate 
Error_Array = xor(Data_Out,Data_In); 
Total_Errors = sum(Error_Array); 
Error_Rate(h,:) = Total_Errors/length(Data_In); % change 1 back 
to h after troubleshooting 


2 


end % ends the outer loop with counter k - for variable noise 


% change Error rate into a horizontal vector - data array 
Error_Rate_PN_RCVR = Error_Rate'; 





sexport error rate 

save ('Error_Rate_PN_RCVR', 'Error_Rate_PN_RCVR') % Saves Error_Rate 
variable into a separate file 

save('SNR', 'SNR') % Saves EbNo_array variable into a separate file 
save ('Eb', 'Eb') % saves the value of Eb - it will be included on the 
plot 





ol? 


Error performance will be plotted in a separate file. 


ol? 


End File 


SUPPORTING FUNCTION ‘UPSAMPLE2’ MATLAB CODE 


ol? 


Name: LCDR Frank Cowan 

Title: upsample2 - a function 

File Name: upsample2.m 

Simulation Project: NPS EE Thesis 

Date Created: 23 Apr 09 
MATLAB/Simulink Version: 7.6.0, R2008a 


AP AP ol? 


ol? 





ol? 


ol? 


This function upsamples a given data sequence by a user-entered 
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factor. 


fe) 


of the 


fo) 


% cyclic PN code for my thesis. 


function [out] = upsample2 (in, U) 
out = []; 
for k = 1:length (in) 


out = [out in(k) *ones(1,U)]; 


end 


% End File 
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% It is designed as a component of PN_Code_RCVR_VF, 


the 


simulation 
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